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Abstract 

Within the framework of the operator product expansion (OPE) and the 
renormalization group equation (RGE), we show that the temperature and 
chemical potential dependence of the zeroth moment of a spectral function 
(SF) is completely determined by the one-loop structure of an asymptotically 
free theory. This exact result constrains the shape of SF's, and implies a 
highly non-trivial functional form for the SF near second order, or weak first 
order, phase transitions. Phenomenological parameterizations of the SF, often 
used in applications such as the analysis of lattice QCD data or QCD sum 
rule calculations at finite temperature and baryon density, must satisfy these 
constraints. 
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Typeset using REVTeX 



Hadronic properties at finite temperature and baryon density are of great importance in 
the phenomenology of heavy ions collisions, star interior and the early universe. Moreover, 
the theoretical expectation of transitions to a chirally symmetric phase and, perhaps, to a 
quark-gluon plasma phase contributes to the interest in studying the effect of matter and 
temperature on the QCD vacuum. 

Our present understanding of QCD at finite temperature (T) and baryon density (or 
chemical potential p) is mainly limited in the Euclidean realm, due to the lack of non- 
perturbative and systematic calculating tools directly in the Minkowski space. Typical 
methods, with QCD Lagrangian as the starting point, are the OPE and lattice simula- 
tions. Because these two formulations are intrinsically Euclidean, only static quantities are 
conveniently studied. In order to gain dynamical informations, which are more accessible 
experimentally, the analytic structure implemented through dispersion relations often have 
to be invoked within the theory of linear response. 

The real-time linear response to an external source coupled to a renormalized current 
J(x) is given by the retarded correlator: 

K(x;T,p) = e(x o )([J(x),J(0)}) T ^, (1) 

where the average is on the grand canonical ensemble with temperature T and chemical 
potential p. Disregarding possible subtraction terms, we can write the following dispersion 
relation for the frequency dependence of the retarded correlation function: 

K{uj,k;T,p) = du — — — — . 2 

Jo u — \ui + ie) z 

For convenience, we discuss only the uniform limit (k = 0). Upon analytic continuation, 
uj —>■ iQ, the dispersion relation in principle connects the correlator in Euclidean region to 
the SF, which embodies all the real-time information. 

Excluding special cases, where the SF p{u; T, p) is experimentally measured, neither 
side of Eq.(^) is completely known. In order to extract physical information in terms of 
elementary excitations, the standard approach consists of calculating the left-hand side of 
Eq.(||) (typical means are current algebras, lattice simulations and the OPE), assuming a 
phenomenological motivated fit to the SF p(u; T, p) and using the dispersion relation, Eq.(^), 
to determine the physically motivated parameters of the fit. This approach has been very 
successful at zero (T,p), due to the fact that we have a good qualitative understanding on 
how to parameterize SF's. Two important example are the QCD sum rules pioneered by 
SVZ [0,0 , and the analysis of lattice QCD data || . 

The QCD sum rules approach has been extended also to systems at finite tempera- 
ture . The lack of experimental data, and of reliable nonperturbative calculations has 
prompted people to use the same kinds of parameterizations that have worked so well at 
zero (T,p) with, at most, perturbative corrections. We believe that physical results at finite 
(T, p) can be strongly biased by this assumption. In fact, recent interpretations of lattice 
simulation data appear to indicate the existence of such problems. 

The purpose of this letter is to derive exact sum rules that constrain the variation of 
SF's with T and p. The derivation, based on the OPE and the RGE, has closer analogies 
to the analysis of deep inelastic lepton scattering experiments than to the QCD sum rule 
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approaches [Q,[|£|-|(| . Moreover, our derivation is not based on expansion for large or small 
(T, /i), making these sum rules valid everywhere in the (T, /z)-plane, and in particular near 
phase transitions. In addition, we apply these sum rules to the chiral phase transition, and 
demonstrate that SF's in some channels are drastically modified compared both to their 
zero T and perturbative shape, at least when the phase transition is second order or weak 
first order, confirming our worries about the non-trivial effect of finite T or baryon density. 

In an asymptotically free theory, the OPE yields, e.g. in the MS scheme, the large-Q 2 
asymptotic expansion 

K(iQ; T, p, k) ~ K (iQ, k)+J2 C n (Q\ g 2 (K), K)([O n ] K ) Tilt , (3) 

n 

where g 2 (n), [O n ] K 's and C n 's are, respectively, the coupling constant, the renormalized com- 
posite operators and their corresponding Wilson coefficients at the subtraction mass scale 
k. It is important to notice that the information of the ensemble average is encoded in 
the matrix elements of the composite operators, while the Wilson coefficients and K are 
independent of T and p. Although the matrix elements ([O n ] K )T,fj, cannot be determined per- 
turbatively, the Q 2 -dependence of Wilson coefficients C n is dictated by the renormalization 
group equation flQfl , and given by 



g 2 {^) 



{l + Otf(Q))}, (4) 



where d n is the canonical dimension of the operator O n minus the dimension of K, and 
c n (g 2 (Q)) is calculable perturbatively. The pure numbers 7« (i = J, n) and b are related to 
the anomalous dimensions of J, O n and to the /3-function as follows 



= - ll9 2 + 0(g*) , and (3 = -bg* + 0(g b ) . (5) 

To study the dependence of K on (T, /i) we only need to consider the difference AK (iQ) = 
K{iQ- T, fx) - K(iQ; T', li') and 



oo 



where Ap(u) = p(u;T,fj) — p(u;T',pf). This subtraction is crucial to remove K (iQ,n), 
which contains all the terms not suppressed by a power of 1/Q 2 , and also to make AK(iQ) in- 
dependent of the renormalization point k. Finite masses give corrections of order m 2 (Q) / Q 2 , 
with m 2 (Q) that runs logarithmically and hence can be ignored, if we are only interested in 
the lowest moment of the subtracted SF. 

At this point we have expressed the left-hand side of Eq. @ as an asymptotic expansion 
of the form: 

AK(iQ) ~ £ ci:){K) ^}} 0nU) [g 2 (Q)r Vn , (7) 

where A([O n ] K ) denotes the difference between the expectation values of [O n ] K in the en- 
sembles specified by (T, p) and (T',p'), while the exponent r\ n and the Q 2 -independent 
coefficients c^(k) are known perturbatively. 
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We proceed by making an analogous asymptotic expansion of Ap{u): 



M-Hfi^SEWwr (8) 

n LI n 



For notational clarity, we have ignored exponentially suppressed terms and the fact that 
there can be more than one rj n and £ n for each n. We then obtain the sum rules by imposing 
that the asymptotic expansion of the right-hand side of Eq. ([]), which we get by inserting 
Eq. (H) in the dispersion integral, matches the left-hand side obtained by the OPE. We 
wish to point out that terms like g 2 T/y/v?, even though allowed by dimensional reasons, 
cannot appear in Eq. (J8|), because they would generate terms like g 2 T '/ ' in the OPE 
series. On the contrary, these kind of terms can be generated by infrared instabilities in 
quantities that do not possess scale separation of the OPE type, e.g. thermodynamical 
quantities. Perturbatively, terms like ln(w/T) would appear in the SF at high-T limit with 
u fixed, but disappear in the large-u limit with T fixed. To carry out the dispersion integral, 
the integration region is splitted into three intervals, (0, A 2 ), (A 2 ,Q 2 ) and (Q 2 ,oo) with 
A 2 large, but otherwise arbitrary. The integration in the first interval can be carried out 
naively, while standard integration by parts techniques have been applied to the last two 
intervals. In a long paper |TT[ we present those technical details that make the procedure 
sketched above rigorous. In particular, we show how to make the asymptotic expansion of 
the dispersion integral (e.g. how Eq. fliTf ) is obtained from Eq. (|^)), and how to match it 
with the asymptotic OPE series of the correlator. At the moment we only need to consider 
the leading terms in the expansions of, respectively, the left-hand side and the right-hand 
side of Eq. (^j) : 



1 

Q 2 



Ap 



M 



-[g 2 (Q)f 



(10) 



In Eq. fllPf), if £ n > 1, Ap can be shown to be equal to the zeroth moment of the subtracted 
SF; note that the zeroth moment of a function whose asymptotic expansion is Eq. (§) is 
infinite, if £ n < 1, since J^dx x -1 (log x)~^ n = oo. 

For the sake of concreteness, let us examine the consequences of matching Eq. (^) and 
Eq. ( |TUD in the case we are interested in, i.e. n = 1 and di = 2. 

If the OPE calculation produces 7] n = 0, then £ n must be an integer greater than one, the 
zeroth moment exists and is given by 



(11) 



du 2 Ap(u) = cf ) (K)A([0 1 ] K ). 
If r) n > 0, then £ n = 1 + rj n , the zeroth moment is again finite and equal to zero: 

roo 

/ du 2 Ap(u) = . 
Jo 



(12) 
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Our main results, Eq. (|TT|) and Eq. (0), can be expressed in physical terms as follows. 
The zeroth moment of a SF for a current J whose OPE expansion yields r] n > is indepen- 
dent of T and //, while the same moment for a current with i] n = changes with T and \i 
proportionally to the corresponding change(s) of the condensate(s) of leading operator(s). 
Furthermore, the zeroth moment clearly is independent of the renormalization point. 

Since £ n < 1 or £ n = 1 in Eq. (|8|) produce terms such as [(7 2 (Q)] _1 or ln[g 2 (Q)] in Eq. (|IT)|). 
this possibility, which implies that the zeroth moment of the SF is infinite, can again be 
deduced from the OPE result (rj n < or the presence of ln[g 2 (Q)]). We remark that, even 
when the moment is not finite, the asymptotic expansion is still well defined and, therefore, 
the whole approach is still correct [ |TT|j . 

At this point several more general comments are appropriate [|l| 



(1) Our derivation relies on the fact that an asymptotically free theory allows a perturbative 
expansion at short distances, making practical the use of the OPE and of the RGE. We 
understand why only short distance physics is involved if we realize that the integral over 
frequencies reduces Eq. (|]) to the ensemble average of the equal-time commutator of the 
currents. Therefore, results such as Eq. (|TT| ) and Eq. (|T2| ) are completely determined by 
the one-loop structure of the theory and the particular current under exam. Within the 
framework of the OPE and of the RGE, the derivation is exact, i.e. there are not corrections 
from higher-loop graphs. 

(2) Flavor, or other non-dynamical quantum numbers, does not change the expansions at 
the one-loop level in an essential way and, therefore, one can derive analogous sum rules by 
using other kind of subtractions, instead of (or in addition to) different T or fi. One such 
example is given by the exact Weinberg sum rules at zero T |13| in the chiral limit. 

(3) The derivation of sum rules for higher moments of the SF requires the complete cance- 
lation of all the lower dimensional operator terms with their complete radiative corrections, 
not just the leading g 2 (Q) corrections; in particular, we also need current quark mass cor- 
rections to the Wilson coefficients. Without appropriate subtractions, higher moments do 
not even converge. 

(4) It is essential to take properly into account the QCD logarithmic corrections when 
deriving the exact sum rules, for the logarithmic corrections not only dictate whether Ap 
satisfy Eq.(|TTl) or Eq.(|i~2"l), but also control the very existence of Ap. This procedure is 
in sharp contrast with the usual QCD sum rule approach, where the convergence issue is 
by-passed by applying the Borel improvement. 

(5) We believe that the (T, /^-dependent part of the leading condensate appearing in Eq. ( |TTD 
does not suffer from the infrared renormalon ambiguity. In fact, only the perturbative 
term K can generate contributions to the leading condensate that are dependent on the 
prescription used to regularize these renormalons. But Kq is independent of T and fi: any 
prescription dependence cancels out when we make the subtraction in Eq. @. On the 
contrary, unless we generalize Eq. (|6|) and make other subtractions, sum rules that involve 
non-leading condensates may be ambiguous. 

(6) In addition to applying our general results to QCD, we have also explicitly verified their 
correctness in a soluble model. A section of Ref. |lTJ is dedicated to illustrate the derivation 
in the Gross-Neveu model in the large-iV limit, where we can either calculate the Wilson 
coefficients, (3- and T-functions in vector and pseudoscalar channels, or obtain the exact SF's 
at arbitrary (T, p), and hence directly verify the sum rules. 
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Now let us specialize to QCD and consider four correlation functions: two involving 
the non-conserved scalar J s = and pseudoscalar J P = ^75^ currents (anomalous di- 
mensions 7j s = 7j p = 1/47T 2 ) and two involving the conserved vector J v = ipj^ip an d 
axial- vector J a = ^J^Jb^ currents (7^ = 7^4 = 0). Since the leading operators (dimen- 
sion four) have nonpositive anomalous dimensions, the two non-conserved current have 
Vn > (27j — 7n)/26 > and Eq. ( [L2]) applies, i.e. the zeroth moments of their SF's are 
independent of T and p. On the other hand, the two conserved currents have r] n = and 



a generalization of Eq. (O) applies [11 1, where the three dimension-four operators with 



zero anomalous dimension appear. Two of these operators are Lorentz scalars: mipip and 
a s G 2 = a s tr(G[ ll ,G tlu ), while the third is the energy-momentum tensor: 

V = -${lviD v + IviD^ - mg^ + 2 (ti{G^G v a ) - -^G 2 ^ . (13) 
The sum rule for the vector (a = 6) and axial-vector (a = —10) currents is: 

du 2 Ap(u) = aA([m#]) + A ^ G ^ + 8A([9 00 ]) . (14) 
in 

These exact sum rules should not be affected by instantons, although the value of the conden- 
sates certainly have instanton contributions. The reason is that the instanton singularities 
in the Borel-plane are located on the positive axis starting at 871" 2 , and, therefore, contribute 
to correlation functions only with higher order terms in 1/Q 2 | |14| |. 

Finally, let us discuss some of the phenomenological consequences of these exact sum 
rules. In the pseudoscalar channel, Ap = implies that, in the broken-chiral-symmetry 
phase, the change of the pion pole induced by T or /i is exactly compensated by a corre- 
sponding change of the continuum part of the SF. This result could in principle be verified 
experimentally, but more practically at present it implies a strong constraint in parameter- 
ization of SF's. For example, the residue of the pion pole can not be taken as a parameter 
independent of the continuum part. Next let us consider the scalar correlation function at 
zero frequency, the chiral susceptibility, 

X (T,fi) = Jd 4 x9(x Q )([Js(x),J s (0)})T,, (15) 

^ 2 Pj^) ; (16) 
u 

which diverges when (T, /z) approaches the phase boundary, provided the chiral restoration 
is a continuous transition. The divergence of the chiral susceptibility near phase transition 
can be only produced in Eq.(^) by singularities very close to the origin, when the exact 
sum rule Ap = is simultaneously taken into account. Thus, when approaching the phase 
boundary in the (T, p)-plane, the threshold of the spectral function vanishes (since there is 
no massless poles in the chirally symmetric phase), and a strong peak develops right above 
the threshold. Because in the chirally symmetric phase the pseudoscalar and scalar channels 
are degenerate, the same would also happen to the pseudoscalar SF. This strong peak in the 
pseudoscalar and scalar SF's, which is intimately connected with the critical phenomenon of 
a diverging susceptibility and correlation length near the phase transition, can be interpreted 
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as some kind of quasi-particle, thus confirming the qualitative picture, originally proposed 



in the context of the Nambu-Jona-Lasinio model ||15|| , of the appearance of soft modes near 
the chiral phase transition. 

We believe that similar result hold in the vector and axial-vector channels, even if our 
conclusions cannot be as strong as in the previous case, because of two main differences. The 
first is that now Ap ^ 0; this problem is not very serious, since we only need that Ap is not 
singular crossing the phase boundary. And this can still be argued by means of the sum rule 
in Eq. (|T4]) and the fact that, for a continuous phase transition, the changes of the thermal 
energy A([6 , o]), and of both condensates A([mijjip]) and A([a s G 2 ]) should behave smoothly 
across the critical line. The most serious difference is that, in these channels, we cannot 
argue on physical grounds that the corresponding susceptibility diverges. Nonetheless, there 
exist lattice simulations [|16| showing that the so-called quark number susceptibility rapidly 
increases in the transition. Thus it is plausible that also in the vector channel the SF has a 
vanishing threshold, and accumulates strength just above it. 

If the chiral restoration turns out not to be a second order phase transition (infinite 
correlation length), but rather a cross-over or weak first order transition (finite but large 
correlation), as the lattice data seem to indicate fll7| , we expect the same qualitative features, 
though less pronounced: a peak develops but it does not actually diverge at the transition. 

In summary, we used OPE and RGE to derive exact sum rules at finite T and p valid for 
asymptotically free theories. We found that, depending on the theory and on the current 
under study, the zeroth moment of a SF is either independent of T and p, or its change is 
related to the corresponding changes of the condensates of operators of lowest dimension. 
The inclusion of the logarithmic corrections in our derivation is crucial, for the logarithmic 
corrections not only determine whether a sum rule should obey Eq . (|TTD or Eq.(|T2"D but also 
control the very convergence of the zeroth moments. These exact sum rules strongly con- 
strain the shape of SF's in every point of the (T, p)-plane, especially near phase transitions. 
We urge whoever parameterizes a SF, e.g. in QCD sum rule type of calculations or to in- 
terpret lattice simulations, to incorporate these exact constraints. For instance, the change 
of the pion residue is not independent of the change of the continuum part. The present 
knowledge about the order of the QCD chiral phase transition together with the information 
carried by the sum rule in the scalar channel has important consequences near the chiral 
restoration phase transition. More specifically, we find that the scalar and pseudoscalar SF's 
differ drastically from their zero (T, p) shape, and develop a strong peak right above thresh- 
old; we argue that the same might be true for the vector and axial-vector SF. Therefore, 
results derived using the same type of parameterization used at zero (T, p) should be taken 
with great caution. In the future, we plan to generalize these results to baryonic currents 
and analyse their phenomenological consequences in greater detail. 

This work is supported in part by funds provided by the U.S. Department of Energy 
(D.O.E.) under contracts #DE-AC02-76ER03069 and #DE-FG06-88ER40427. 
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